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Abstract 

We renormalize QCD to one loop in coordinate space using constrained differential 
renormalization, and show explicitly that the Slavnov-Taylor identities are preserved 
by this method. 


Differential regularization and renormalization Q is a method that works directly 
on Feynman graphs in coordinate space, substituting singular expressions by derivatives 
of well-behaved distributions. It has proved to be quite simple and convenient in a 
number of applications. In gauge theories, however, a problematic feature arises: the 
Ward identities must be studied explicitly to fix the arbitrariness of the method in such 
a way that gauge invariance is preserved. 

A solution at the one-loop level is the symmetric procedure of differential renor¬ 
malization proposed in Ref. |||. This so-called constrained differential renormalization 
(CDR) introduces no ambiguities and has been shown to preserve Abelian gauge in¬ 
variance g! and supersymmetry in a non-trivial calculation in supergravity 0j. It 
is the purpose of this letter to apply CDR to a non-Abelian gauge theory and study 
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the corresponding Slavnov-Taylor identities We shall consider a Yang-Mills theory 
with gauge group SU(N C ) coupled to Nf Dirac fermions in the fundamental repre¬ 
sentation, i.e., QCD with N c colours and Nf quark flavours. The calculation of the 
gluon selfenergy and the triple gluon vertex, using conventional differential renormal¬ 
ization, was carried out in Refs. [i|] and Q, respectively. The background field method 
|?[] was employed there because it allows a much more direct determination of the (3 
function and leads to simpler Ward identities. We use the conventional formalism in¬ 
stead, precisely for the last reason: we want to test CDR in the most complex case, 
and show that it preserves the (more involved) Slavnov-Taylor identities rather than 
the Ward-Takahashi like identities of the background field formalism. Nevertheless, 
for comparison with Refs. 01, we have also applied CDR to the background field 
calculations mentioned above. 

In Ref. H it was argued that CDR preserves the Ward identities of Abelian gauge 
invariance because it maintains the properties that are required for their derivation, like 
the fulfilment of equations of motion for renormalized expressions or the commutativity 
of differentiation with renormalization. Actually, this argument applies equally well 
to the case of non-Abelian gauge invariance, since the structure of the interaction 
Lagrangian is never used. Of course, the symmetry of the Lagrangian is essential for the 
fulfilment of the corresponding Ward identities, but this is a matter of combinatorics 
Pj and CDR does not interfere with it (essentially, CDR ensures that the building 
blocks of such combinatorics behave correctly). A new feature of the non-Abelian 
case is the appearance of composite operators in the Slavnov-Taylor identities (in the 
Zinn-Justin form |J, which we shall use here). In Ref. [jlj] it was shown that differential 
renormalization can be directly applied to diagrams with operator insertions. The same 
holds for the constrained procedure. 

In the following, after writing the Lagrangian, we give the renormalized expressions 
in coordinate space of all the singular one-loop 1PI Green functions of elementary fields. 
The diagrams with operator insertions that contribute to the Slavnov-Taylor identities 
have also been calculated, but the explicit results are not given here. Then, we write all 
the Slavnov-Taylor identities involving these renormalized Green functions. We have 
used a symbolic computer program to verify that they are indeed fulfilled. The QCD 
Lagrangian in the Feynman gauge, written in Euclidean space and including ghost 
terms, reads 

1 1 Nf 

£ = + d^iD^r + £ HP + , ( 1 ) 

2=1 

with 

= d,A a v - d v Al + gf abc A b ^Al , (2) 

{D^) a = d^r, a +gr bc Ay, ( 3 ) 

D^ l = <VR - i gA^T^i . (4) 

T a are the SU ( N c ) generators in the fundamental representation, and f abc the structure 
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constants. This Lagrangian is invariant under the BRST transformation [ 101 


SAl = 

-(D^rs A = 

sAl5X, 

fj, 

(5) 

5fj a = 

-d^dX = srj a 5X , 

(6) 

5y a = 

~f abc V b rj c SX 

= ST] a dX , 

(7) 

SVi = 

-gr'tfViSX e 

= sTjdA, 

(8) 

S^i = 

-g*iT a ri a 6X = 

e s$i6X , 

(9) 


where 5X is a constant Grassman parameter. In order to obtain simple Slavnov-Taylor 
identities for the 1PI Green functions, it is convenient to add to the Lagrangian source 
terms not only for elementary fields but also for their BRST variations Q: 

/ _ Nf 

Sources = - I + fj a £ a + + £(*iXi + Xi*i) 

N f 

+K;(sA “) + L a (sr , a ) + 

2—1 

The sources of BRST variations are mere spectators in the Legendre transform defining 
the generator functional of 1PI Green functions. The coordinate space Feynman rules, 
including those for insertions of BRST variations, are displayed in Fig. [I]. 

The singular 1PI Green functions of QCD at one loop are the gluon, quark and 
ghost selfenergies, and the quark, ghost, triple gluon and quartic gluon vertices. The 
procedure to calculate them with CDR is simple j|, [ij. First, one expresses each 
contributing diagram in terms of the basic functions defined in Ref. j3j]: 



A m = A m (x)8(x ), (11) 

B mi m 2 [0\{x) = A mi (x)O x A m3 (x ), (12) 

Tmim2?n3 [(D\(x,y) = A mi (x)A m . 2 {y)O x A m3 (x + y), (13) 

Qmi)7i2m3rri4 [^] (^h Vi ^0 = A mi (x) A m2 (y) A m3 (z)0 A m4 (x T y T z) , (14:) 


where A m (x) = An - g ^j ie F e y nman propagator in (Euclidean) coordinate 

space, with K\ a modified Bessel function, and O is a differential operator. We shall 
suppress mass subindices when all of them are zero. Then, the singular basic functions 
are replaced by their renormalized expressions, which can be found in Tables 2, 3 and 
4 of Ref. ||. For brevity, we refer to those tables and do not reproduce here the 
renormalized basic functions that appear in our calculations. Also, since the procedure 
of CDR has been explained in detail in Refs, [pj, |3|, we shall directly give the final result 
for each Green function. The Feynman diagrams contributing to 1PI Green functions 
of elementary fields can be found, e.g., in Ref. [11]. 


The renormalized gluon selfenergy in terms of renormalized basic functions reads 


<A“(S!)4(X 2 ) >= g 2 S ab 


N r 


(3d,,dv 


25^a)B R [l] 


4B R [d,A] 
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(15) 


N f 


+ 


Y. [( — 29 m 9j, + 5^(0 + 2m-)j B^. m .[l] - 2B^. m .[n] + 4B^. m . [9^9,, 


2—1 

where all functions, here and in the rest of two-point functions, depend on the coor¬ 
dinate difference x = x\ — x 2 and, as prescribed by differential renormalization, total 
derivatives are supposed to act formally by parts on test functions [j|]. Throughout this 
paper, < <f>i.. .<t> n > represents the one-loop correction to the 1PI Green function of 
the fields (f>\... (j) n . That the expression in Eq. (^) is transverse, as required by gauge 
invariance, is only apparent when the renormalized basic functions are substituted by 
their explicit expressions in the mentioned tables. In the case of massless quarks the 
result is quite compact: 


< K (®i)^(* 2 ) >= - 


1 


gH ab (9^ - S^n) 


1447T 2 


1 lo {tx 2 M 2 

(151V C - 61V/) —D -^-2 -+ (2N C - 2 N f )5(x) 


(16) 


M is the renormalization scale. The dependence on the renormalization scale of Eq. (15) 
is 


M 


9 

9M 


1 


< A“( Xl )K(x 2 ) >= — g '5 ab (5N c - 2 N f )(d»d v - 5^a)6(x). (17) 

In the background field method, the background gluon selfenergy renormalized with 
CDR reads (for massless quarks) 

< B a J Xl )Bt(x 2 ) >= - 1 


144^*' 

/OOAT /.jy \ 1 10g.T 2 M 2 

(33 N c -6N f )^n - — 


+ (2N C - 2N f )5(x) 


(18) 


and differs from the result (for Nf = 0) in Ref. Q by a finite local term (which 
can be absorbed into a redefinition of M). The one-loop j3 function of QCD, j3 = 

3 _ 

— ^fp-(lliV c — 2 Nf), can be directly read from the scale dependent part of Eq. (18). 

The quark selfenergy is proportional to the corresponding Green function of 
QED. In terms of basic functions, the renormalized quark selfenergy reads 


< tff(si)#f(s 2 ) > = g z S 


2 X AB N c ~ 1 


R 


N r 


1 g 2 5 ABN <~ 1 


647T 4 


~ 2 m,;B 0m 

(0 — 4m 


■W) 


N, 


, n 2 J < o{m i x) 
(□ -mj- 


M 1 


+ 2tt log — j5(x 


mt 


+ m 2 K 0 (rriiX ) 

x z 


(19) 


where Kq is a modified Bessel function [^] and M = 2M/'je, with 7 # = 1.781... the 
Euler’s constant. The scale dependence is 

9 - A , x ^ b t - 1 


M 


9M 


< (x 2 ) > = 


167T 2 


]\J 2 _ 1 

g 2 § AB _c_ -^ . 


( 20 ) 
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The renormalized ghost selfenergy is 


< V a (xi)ff(x 2 ) > = — -g z N c 5 ab OB u [1] 


1 log X 2 M 2 

128 vr 4y x 2 


and its scale dependence, 


M WM < = ~-^23 2 5 abN c^{x) 


( 21 ) 

( 22 ) 


For the three-point functions we use the shifted variables x = x\ — x 3 , y = x 3 — x 2 - 
Unless otherwise specified, all basic functions are assumed to depend on these two 
variables in the following. To avoid too lengthy expressions, we shall only give the 
final expressions in terms of basic functions. The renormalized expression of the quark 

vertex is 

< ^(sr)*f (z 2 )A“(x 3 ) >= -i g\T a ) BA ji- [(—2m*9+ + m ? 7 m 

+ fi Xr hxfi V ) TmimfO[l] + (4mi<5 Ma - 7a7/i^ - T miTOi o[9 a ] 

- 7mT* mi0 [°] + 2T* .^oM] + iV c gm*(^ 7/i + 7^)^ [1] 

+ Q(0 x 7/z7a + lain ft 1 ) + 7/^ - 2 7a 9+ - <V*(3m* + 2^+)j T 00mi [9 a ] 

+ 7/uToomj [ n ] + 2T^ m . [9 M | , (23) 

where we have introduced the notation 9 + = 9"+9 2/ . The scale dependent part reduces 
to 

M m < Ufe)*?feMJN) >= ^/(r 0 )®- 4 (3AT C - T) 7 ^(i)<5(!/) ■ (24) 


For the ghost vertex we have: 

< ? ? 6 (.ti) ? 7 c (.t 2 )A“ (x 3 ) >= \g*N c r bc {d* • 9^T[1] 

+ («□* - 29" • 9*) + 9;9£ + 9^9-) T[9«] + 9^T R [D]} , (25) 

where 9“ = d x — 8 y and d x ■ d y = d x d y . The scale dependence is 

M ^M < r ? ft ( Xl )^ C ( ;c 2 )A“(x 3 ) >= j^ N cf abc 5(xi ~ x 3 )d X2 5(x 2 - x 3 ). (26) 

Here, we have come back to the original variables, to make explicit that it has the same 
form as the corresponding term in the Lagrangian. 


5 



The renormalized expression of the triple gluon vertex is quite large, even in 
terms of basic functions. We split it into the pure gauge part and the fermionic part. 
The results are 


< > G = g 3 f abc N c 

^S up d~ (B[l](x)<5(x + y)) + ^{Sp V d x p - 5p P %)B[l](x)S{y) 


+1 ( 3 %%d~ + s^{u x + 3d x • a»)% + s pp (5n x % - 2 ( 0 y + 2d x ■ &>)%)) t[i] 

+ \ (<W^ - 3d x d x ) + 5 fJ-p (8d x d x - \&idl - 3dldl) 

+ <5 pQ (8df,d x - 2 Offfi - 3d x d x + 5p U (a x - 7d x • &>)) 

+ Spa (6%% - 13%% - dl&l - 3d x d x + 8 vp (U* - 10D 1 ' + 10^ • &•))) T[d a ] 
+ \ (2 Sp V d x + 5p P (2% - 5d x )) T*[D] + Ad x T R [dpd u \ + 4c^T R [d M cy 


+ 5p U d x T H [d a d p ] + SppdiT H [d a d u \ - 4T K [dpd u dp] 

+ 

for the pure gauge and 


x y 
H <-> v 


(27) 


N f 


< 


Al( Xl )AUx 2 )A c p ( X3 ) > F = 2g 3 f abc J2 [ml (s, w d x - 6 pp d~) T mimimi [l] 

i =1 l 

+ [Sp V d x p dl - 5 pp {d x v dl - %%) + l -5 pa {%% - %% - Sp V {d x • &> + mf)) 

+ Spa (-%% - %% + s up (d x • d“ - rn?))) T m imi mi [d a ] 

- (Sp U dp - Sp P d~ ) T* mimi p\ + 2 d x p T R imimi [dpd u \ + 2d^T R . m . m .[d p d p ] 


2Sppd a T? mirnirni \d a dv\ + r,Spv'Tm-im-imi IP^p] A Sp P T mim . m .[Dd u \ 


1 


R 


t R 


- 2T 


R 

mifnim■ 


,[dpd u d p 


+ 


2 

x y 

fi u 


(28) 


for the fermionic part. The Bose symmetry of the three gluons is not obvious because 
of the use of the shifted variables. Only the symmetry under interchange of A a p (xi) and 
A b u (x 2 ) (corresponding to ( x,y,a,g ) (—y,—x,b,i/)) is explicit. The complete scale 
dependence, written in the original variables, reads 


M 


d 


— < A“( Xl )Al(x2)A C p(x 3 ) >= —yr^N c - N f ) 

- d?) + M^ 2 - d x n + Spp(d x 3 - dA [S( X1 - x 3 )S(x 2 - a*)] .(29) 
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In Ref. the bare triple gluon vertex in the background field formalism (with 
massless quarks) was found to be conformal invariant for non-coincident points, if 
the Feynman gauge is employed. This fact, together with the fulfilment of the Ward 
identities, implies that this three-point function must be a linear combination of the 
two permutation odd conformal tensors D s ^ m {xi,x 2 ,x 3 ) and C'^p im (x i,x 2 ,x 3 ) |l3| . 
The explicit expression in terms of these conformal tensors was also found in Ref. Q. 
Differential renormalization was then used to treat the singularities at coincident points. 
Of course, renormalization breaks conformal invariance, but the Ward identity relating 
the triple background gluon vertex to the background gluon selfenergy was enforced by 
adequately adjusting the renormalization scales that appear in the process. We have 
applied CDR to the renormalization of the conformal tensors D symm and (7 s y mm ; and 
checked that the resulting amplitude directly fulfils the Ward identity (if the CDR result 
for the gluon selfenergy is used) 1 . No adjustment is needed a posteriori. The tensor 
£<symm but ambiguous, and the result in CDR differs from the one given in 

Ref. [jo] by a finite local term. In our case, it gives a non-vanishing contribution to the 
Ward identity. The result for D symm also has an extra local term with respect to the 
final one in Ref. [6]. The discrepancies are due to the fact that, while in CDR everything 
is fixed from the start, in conventional differential renormalization the renormalization 
scales can be adjusted, in general, in more than one manner to preserve the Ward 
identities. 


Finally, we have calculated the quartic gluon vertex, but the final expression 
is too lengthy, even in terms of basic functions, to be written here. We only give the 
scale dependent part: 


M 


d 

DM < ^ 

1 


< A a Jxi)A b l/ (x 2 )A c Jx 3 )A d a {x 4 ) >= 


abr rcdr / 


2l^ 9 {Nc + 2Nf) [f f V*’ 5 '" ~ <W<W) 

+ f acr f dbr (5^ p - VW + f adr f bcr {^^ P a - S pp 8 ua ) 


(30) 


The renormalized 1PI Green functions we have calculated satisfy the renormaliza¬ 
tion group equation 

d d d Nf 

M dM + + lmiTni d^- ~ nA1A ~ n ~ ^ n * i7 * 1 

y 1 i =i 

where ua, n v and nq, i are the number of gauge, ghost and Tquark fields, respectively. 
The coefficients can be easily obtained from the scale dependent parts given above; the 
standard values are recovered: 

_ P= - - 2N f ) > ( 32 ) 

Tt is important that we have used the form of the conformal tensors given by Eqs. (4.1) and (4.16) 
of Ref. [!6). The triple gluon vertex can be expressed in terms of the conformal tensors in this form 
using just the Feynman rules and the Leibniz rule for derivatives, which is an allowed operation in 
CDR. Hence, no ambiguity is introduced in the process previous to our calculation. 


p(n A ,n^,ri5, 1 ,...) _ Q ^ ^) 
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7A = 


(33) 

7n = 

- 92 JV 

32vr 2 c ’ 

(34) 

7 = 

9 2 N*-l 

(35) 

32vr 2 N c ’ 

7 mi = 

3 2 iV c 2 - 1 

(36) 

I 67 r 2 ^ N c 


Note that N c and are the SU ( N ) Casimir invariants of the adjoint and fundamen¬ 
tal representation, respectively. Unlike in the background field method, the anomalous 
dimension of the gauge field is not directly related to the /3 function. The same j3 
function is obtained from all vertex functions, showing that there is a single coupling 
g. This is a consequence of the Slavnov-Taylor identities for the scale dependent parts 
of the Green functions. 


Let us now write the full set of Slavnov-Taylor identities for the complete 1PI Green 
functions. They are a bit more involved than the ones for connected Green funtions 
[14|. The general form in terms of the effective action can be directly derived from the 
BRST symmetry of the Lagrangian |J. Using the source terms given by Eq. (0) and 
suppressing all indices it reads 


/ 


d 4 z 


'sr sr 

JaJk 


st sr sr sr 
JqlL + JnM 


TT ST ST „ 1 

Tt”TaT ^ 

h\k SN Sg 


= 0 . 


To one loop, the quadratic terms can be “linearized”: 


(37) 


ST ST _ (JT(°) <5rW (STW <5T(°) 
S<f> SJ sc f, S4> 5J S <!> + S(j) SJ s< j, ' 

On the other hand, the ghost equation of motion, 


n 5T , ST _ A 
^ SK% Sfj a ’ 


(38) 


(39) 


allows to simplify the identity for the gluon selfenergy. This equation is trivially ful¬ 
filled to all orders in any renormalization scheme commuting with differentiation and 
preserving the structure of the Lagrangian. Writing explicitly the tree-level pieces, the 
one-loop Slavnov-Taylor identities read (we name each identity after the Green function 
with the largest number of elementary fields) 


• gluon selfenergy identity: 


o = d; < A ;(x)A b u (y) > , 


(40) 
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quark vertex identity: 


0 = 91,(T a ) AB 5{x -y)< y b (z)K-(x) > (x^f^A^z) > 

- {@ y + rm) < Tf (x)y b (z)N A (y) > - < rj b {z)^f(y)Nf(x) > (fl x -m*) 

- g(T b ) CB 5(x -z)< (x)^ A (y) > 

+ g(T b ) AC 5(y - z) < (s)ttf (y) > , (41) 


y/iosi vertex identity: 


0 = 


d x L < A^(x) V c (y)fj d (z) > +gf acd d ;; (<5(y - z) < r/>(x)A» > 

+ D z < y b {x)y c {y)L d {z) > +gf abc 5(x - y) < y a (x)fj d (z) > ,(42) 


b <-> c 
x <->• y 


triple gluon vertex identity: 


0 = 


y/ adb <% - z) < A c v {x)A a (y) > -□* < r, b (z)A d (y)K c u (x) > 


< rf(z)ff(x)A d (y) > 




c <-> d 

> 

+ 

v <-> p 



X <-> y 


+ gf acd (-6^,(2% + d y ) + 6 vp (d£ - d y ) + 6 pp (d x + 2d y )) 
(<5(x - y) < y b (z)K d (x) >) - d* < A b ^z)A c v (x)A d (y) > , 


(43) 


quartic gluon vertex identity: 


0 = 


9f acd (S„(2d y + d*) + 5 up (d y - d z p ) + 5 pp (d y + 2dl)) 

(5{y -z)< r, b (t)A*(x)K“(y) >) + gf bae 5(x - t) < Al(y)A d (z)A a a (t) > 
~ □* < Al{y)A d {z)y b {t)K{x) > -<% < Al{y)A d Jz)y b {t)f] e {x) > 



e c 


e«(i 

+ 

G V 

+ 

er <-» p 


x y 


x z 


+ 9 2 ( f acr f der (S pp S U(T - 5 pa 5 up ) + f adr f ecr (5 pa 6 up - 5 pu 5 pa ) 

+ f aer f cdr {5 pv 5 pa - VM) (<5(* - V)5(x -z)< y b {t)K d {x) >) 
- Si < A b M)AUy)A d Jz)A e a (x) > . 


(44) 


The diagrams contributing to Green functions with insertions of BRST variations are 
depicted in Fig. |2[ Verifying that the renormalized Green functions calculated above 
satisfy these identities is not straightforward, due to the length of the expressions 
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involved and to the fact that we are not dealing with simple functions but with distri¬ 
butions. The obvious way to compare distributions is to make them act on a general 
test function. In particular one can perform a Fourier transform without loss of in¬ 
formation^. This has the advantage that the resulting finite integrals can be treated 
with standard momentum space techniques and that, for non-exceptional momenta, 
no infrared divergencies appear. We have used a Mathematica -based program to carry 
out the algebraic operations, linked to LoopTools [|l5|] , which calculates numerically 
one-loop integrals. As advanced, CDR respects the Slavnov-Taylor identities. 


Summarizing, we have applied CDR to the one-loop singular 1PI Green functions 
of QCD and have verified that the Slavnov-Taylor identities are preserved. The exten¬ 
sion to Yang-Mills theories with a more general gauge group does not introduce new 
complications, as far as renormalization is concerned, so CDR should treat the general 
case equally well. Here we have sticked to SU(N c ) because the computer implementa¬ 
tion is simpler [|l6|. One could also wonder about the performance of CDR when the 
gauge symmetry is spontanously broken. Again, as CDR does not depend on the struc¬ 
ture of the interaction Lagrangian, there should be no extra problems in dealing with 
this case. The authors of Ref. [|D|, using the programs described there, have checked 
that CDR renders a transverse vacuum polarization in the electroweak standard model, 
and recovers the standard physical results for Z-Z and W-W elastic scattering. These 
examples also test the inclusion of scalar fields in a non-Abelian gauge theory. Finally, 
let us point out that very recently it has been found (in momentum space) that CDR 
and regularization by dimensional reduction produce equivalent results at the one loop 
level |15| . This gives an alternative explanation for the preservation of gauge invariance 
in CDR. 
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b 


[1 a b V 

x vAA/VWW y 


5 flu 5 ab A(x - y) 


a 


x.>■ 


S ab A(x — y ) 


y (@ x + m)A m (x-y) 




- g 2 [f abr f cdr ^ P S ua - S p(7 S up ) 
+ f acr f dbr (5 P a5 up - 
+ f adr f bcr (S P J P a ~ <WM] 





/ V. 


-gf abc 


il 


c 


il 


c 




-g(T a ) BA 


Figure 1: Coordinate space Feynman rules for QCD, including insertions of BRST 
variations. In vertices, the derivatives (with respect to the vertex space-time point) act 
on the field indicated by the superscript. 
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Figure 2: One-loop Feynman diagrams contributing to 1PI Green functions with in¬ 
sertions of BRST variations. The Feynman diagrams for singular 1PI Green functions 
of elementary fields can be found in Ref. jll|. The diagrams contributing to the (fi¬ 
nite) gluon-gluon-ghost-antighost function, which appears in the quartic-gluon-vertex 
identity, are identical to diagrams fl-f-4, but changing K by an antighost external leg. 
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